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X, INTRODUCTION 

The launch of L&NDSAT has for the first time provided the water re- 
source manager and practical hydrologist with broad prospects for effi- 
cient acquisition of essentially real-time data. These are usable for 
hydrologic land use assessment, surface water inventory, and for the ex- 
traction of information pertinent to soil properties. This information 
has value not only by and In itself, but also to construct the water- 
shed transfer function for hydrologic planning model's, aimed at estima- 
ting peak outflow from rainfall inputs, 

The reduction of satellite data to practical, operational informa- 
tion requires a clear, easily applicable methodology for converting 
these data into quantitative hydrologic parameters. 

The fundamental objective of this effort is the development of 
such a methodology and its transfer to hydrologic users. It was rea- 
lized that such technology transfer could be made far more effective 
by the parallel development and eventual demonstration of the results 
of a model, specifically structured to take full advantage of the ca- 
pability of LRNDSAT - for example, its frequent recurrence and conse- 
quent ability to determine seasonal variations in the watershed's con- 
ditions, The category of planning models was chosen for development 
and demonstration because of its great practical importance in the de- 
sign of waterworks, because of the wide diffusion of such models down 
to capillary levels within the hierarchy of water resources users, 
and because their implementation is relatively simpler than 
■management models, thus making optimum use of the resources available 
for this effort. 
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Consequently, the effort was structured along two major routes: 
the development of a hydrologic planning model specifically based upon 
remotely sensed inputs, including its test and verification from ex- 
isting records; and the application of . LANDSAT data to supplying the 
model's quantitative parameters and coefficients. Included was the 
investigation of the use of LANDS AT data as information inputs to all 
categories of hydrologic models requiring quantitative surface para- 
meters for their effective functioning. 

The effort thus far has consisted of three phases, The first fo- 
cused on the definition of the "drivers 1 ' - those hydrologic processes 
to which peak runoff is most sensitive - and upon the synthesis of a 
simple yet effective model for the estimation of long-recurrence out- 
flows, The results of the first phase effort were presented in the 
Final Report, "The Application of Remote Sensing to the Development and 
Formulation of Hydrologic planning Models," dated January, 1975, The 
second phase extended this work to include the development of a routing 
model for use in sensitivity analyses , and a quantitative investigation 
of the accuracy and completeness of the hydrologic information which 
can be extracted from remotely-sensed Imagery. 

These findings were reported in the January, 1976, final report. 
The current phase has concentrated on validation of the model and 
upon the synthesis of a simple methodology for performing hydrologic 
analyses from. LANDS AT imagery. 


IX. SUMMARY OF RESULTS AND CONCLUSIONS 

The objectives of the current effort have been to couplets and test 
a hydrologic planning model making maximum use of remotely- sensed inputs? 
to document procedures for extraction of pertinent parameters from satel- 
lite data? and to expose the model and techniques to prospective users. 

In the course of this investigation* several findings and conclusions 
have been reached and are described below t 

The Importance of Surface Parameters to Hydrologic Models 

Its surface-related parameters used by existing hydrologic models 
(Manning's S.C-S. ' Curve Number, etc.) can be derived from. LANDS AT 

imagery to the required accuracy. The opportunity Is available to take 
advantage of the radiometric and multi-temporal properties of the data 
to facilitate identification and classification. Moreover, potentials 
exist for remote monitoring of soil and soil moisture which would be of 
great value to hydrologists. 

Planning Model Development 

A simple model which is amenable to remote sensing inputs has been 
developed and tested; It exploits LANDSAT data for surface cover input 
and can accept soil moisture information as this becomes possible. The 
test was performed for the Bucklodge Branch, Maryland, watershed. Ac- 
curacy achieved . to date appears commensurate with those of conventional 
models- The principal purpose envisioned for the model and benchmark 
test is to demonstrate to prospective users that satisfactory results 
can be achieved and to show how similar analysis procedures might be 
applied to their watershed with their models. 


Synthesis of Analytic Procedures 


Our experience has shown that a significant share of hydrologic 
modeling is being performed under the auspices of local governments. 
Consequently, two impediments to the use of satellite data arise. 

First, there is a lack of standard procedures for extraction of hydro- 
logic data. Second, the expense involved with current computer/ optical 
practices is a deterrent to use. This report proposes a procedure for 
simple visual analysis of LANDSAT images which is accurate enough for 
use in conjunction with hydrologic models. It has been determined 
further that these techniques extend beyond the field of hydrology and 
can have equal applicability to agriculture, land use, etc., analyses. 

Contacts with Prospective Users 

To date, contacts have been made with both foreign and domestic 
users. Due to the relatively new state of modeling overseas, the 
response has been greater there. One Italian region has embarked on the 
analysis of a major watershed, spawned in part by this investigation. 
Other regions have expressed interest. Domestic acceptance will be 
fostered by the demonstration of concise procedures and satisfactory 
results. Such a purpose is intended for this document. 
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III. REVIEW OF PREVIOUS 5HASES 

Of critical concern to water resources, planners and engineers is 
the ability to forecast peak flow events. The. capacity to estimate the 
■magnitude and duration of large-recurrence outflows has a significant 
impact upon the accuracy of sizing and designing waterworks, and thus on 
their cost. 

The tool available to the planner for these purposes is the hydro- 
logic model. Although, the inputs of different models vary, all require 
significant quantities of physiographic and hydrologic information? these 
data are typically expensive to obtain and are often only partially avail- 
able, Remote sensing offers a new source of information which formerly 
had to be acquired by less efficient means or ignored altogether, 

The first phases of this effort conducted from February, 1974 to 
December, 1975, addressed four pertinent topics t 

11 Identification of the ’'drivers' 1 of peak flow events, i.e. , 

. the hydrologic phenomena (infiltration , antecedent soil, 
moisture, etc.) to which the watershed's outflow is most 
sensitive, 

2) The development of a model compatible to the maximum degree 
with remotely-sensed hydrologic inputs, 

3) Preliminary verification of the model for actual watersheds, 

41 Initial, identification of the efficiency of remote sensing , 
in determining the parameters of the model. 
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3*1 Investigation of Driver Phenomena 

The purpose of this investigation was to achieve a valid statis- 
tioal comparison of the rates and magnitudes of the hydrologic pro- 
cesses contributing to the runoff from long-recurrence events. This 
was accomplished over a significant sample of watersheds, with wide 
variation of climatology, terrain , and physiography. The comparison 
allowed the determination of which are important and which can be neg- 
lected without significant loss of accuracy, 

Kain falling on a watershed is subject to several processes which 
abstract water and govern the flow. Those which produce the most sig- 
nificant changes to flow rates and volumes '‘drive” the basin outflow. 
Table 1 describes the hydrologic processes? Table 2 presents the 
’'drivers" of each. It shows that several processes can be omitted in 
the formulation of a peak rate model because of their limited inpact. 
For example, the rates corresponding to interflow, percolation, and 
evapotransplration are very slow in comparison with other processes 
such as rainfall, infiltration; and overland flow. Also, interception 
and depression storage become saturated early in large rainfall events 
and therefore are inconsequential to peak event modelling, 

Therefore, except for very special circumstances, peak flow can be 
adequately modelled by considering only precipitation., infiltration and 
surface flow - both overland and in the channels. The means by which 
this can be accomplished will be discussed in Section V. 


TABLE 1 


SUMMARY DESCRIPTIONS OF HYDROLOGIC PROCESSES 


HYDROLOGIC PROCESS 


DESCRIPTION 


Interception 


Infiltration 


Moisture caught and stored on plant leaves and 
steins or other impermeable objects; eventually 
evaporated back into the atmosphere. 


Downward movement of water from the surface into 
the soil. 


A) Interflow 


Lateral subsurface water movement toward stream 
channels. 


B) Percolation. 


Downward movement of water through soil to 
groundwater (area where. pores of soil or rock 
are filled with water) . 


C) Base Runoff 


Water from interflow and percolation which moves 
underground to the channel. 


Evapotranspiration 


A) . Evaporation . 


B) Transpiration 


Precipitation Excess 


A) Depression 
Storage 


B) Surface Flow 


Upward movement of water in gaseous state from 
the surface. 

Movement of water through plants to the 
atmosphere, 


Retention of excess rainfall in surface depres- 
sions. 

Uninfiltrated water which flows' over land 
surface to stream channels. 


C) Channel Flow 


Total Runoff 


Flow of water in natural channels. 

Sum of runoff from underground processes (base 
runoff) and overland flow (direct runoff) . 























TABLE 2 

potentially bsfortant drivers as related to hydrologic processes 


HYDROLOGIC PROCESS 

PRINCIPAL DRIVERS 

Overland Plow 

Slope 

Roughness of Soil 
& Cover 

Drainage Density & 
Pattern 

Infiltration 

Soil Permeability 
Antecedent Soil 
Moisture 

Soil Moisture Capa- 
city 

A) Interflow 

Soil Permeability 
Subsurface Moisture 
Gradient 

Plow Length* Slope 

B) Percolation 

Soil Permeability 
Subsurface Moisture 
Gradient 
Soil Depth 

1 1 ”1 
. Evapotranspiration ' 

Temperature 
Antecedent Soil 
.• Moisture 
Soil Permeability 

A) Evaporation 

j 

1 

B) Transpiration 

Temperature 
Solar Radiation 
Vegetative Cover 
Antecedent Soil 
Moisture 1 

Depression Storage & 
Detention 

Depression Density 
Cover Retention 

Interception 

Duration of Rainfall 
Intensity of Rain- 
fall 

Cover Composition* 
Age* Density 


SECONDARY DRIVERS 



Vegetative Cover 
Slope 

Water Turbidity 
Temperature 



Water Turbidity 
Wind 
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3.2 Development and Validation of Remote Sensing Model 

The following criteria were followed in the structuring of a peak- 
rate model i 

1) The model would be modular, to allow the user flexibility of 
application, 

2) The model would maximize potential input from remote sensing. 


During the first phase a model of the Rational Formula type was 
developed. It is of the formt 

1 


Q - 2L1 5 


£ Cl) (hi 


.6 


5 ** S ,3 C3600) J 


0.4 - — 
*3 


(1) 


Where » 


3 2 

Q « maximum outflow rate , m /sec/km 


L - channel length, m 
1 « length of average strip, m 
n * average banning* s H n'’ 

<(> m routing factor 

K - K T a 2 

K - infiltration and spatial correction factor 
T * rain recurrence interval, years 
s * average slope, m/m 
a l' a 2 r a 3** con * tant *f function of location 


The output of this model is the expected peak discharge rate for a 
watershed for a rainfall of T years recurrence interval (A complete 
derivation is contained in the January, 1975 report)* This model was 
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verified against records of 31 experimental watersheds of the Agri- 
cultural Research Service, The results shewed this model to equal 
or better the performance of other formulations of the Rational type. 

In the second phase* title two criteria of modularity and remote 
sensing input maximization were applied to the development of a more 
comprehensive model capable of routing water overland and through the 
channels. This model as envisioned would have maximal utility for un- 
gaged f predominantly rural watersheds* i»e, , those where remote sensing 
could play a significant role. Completion and verification of the model 
occurred during the current effort and will be reported later in this 
document, 

3.3 Identification of the Efficiency of Remote Sensing in Determining 
Parameters of Hydrologic Models 

Thi3 task was comprised of hydrologic land use ..analysis of water- 
sheds located as given below s 
11 Cftickasha, Oklahoma 
2} oxford, Mississippi 
3} Blacksburg, Virginia 
4J Muddy Branch, Maryland 

The findings of this analysis are summarized in the following! 

11 Substantial hydrologic information can be measured from 

LANDS AT imagery, The parameters identified and measured weret 
surface water bodies, surface cover classes, channel lengths, 
and watershed area. Both single-band and color composite im- 
agery were examined; the latter offered much more information 
and greatly facilitated analysis , 
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2J Of particular value were seasonal. LANDSAT observations . For 
exanple, the October imagery showed the watersheds when vege- 
tation density was low? this made the higher-order streams 
visible. Wore channels could be measured from the Chickasha, 
Oklahoma , image than appeared on the U.S.G.S, topographic map 
at the same scale, LANDSAT imagery clearly shows that ef- 
fective drainage density varies cyclically throughout the year, 

31. It became apparent that published records do not reflect cur- 
rent surface cover conditions, since the watershed^ land use 
typically alters with time, The October, 1973, LANDSAT imagery 
demonstrated that surface cover for the Chickasha Basin had 
changed front that described in the latest published data (1967) * 

4L The Oxford and Blacksburg watersheds were used to evaluate the 
discrimlnability of stream courses, forest cover, surface water 

t 

and impermeable areas from LANDS AT images, The relative value 
of each of the four MSS bands was assessed; Sand 7 was found 
moat useful for surface water; while Band 5 offered more detail 
' in vegetated areas. It was concluded, though, that composite 
imagery optimized the contributions of all bands. 

51 A complete hydrologic land use analysis was performed using 

color composite imagery for the Muddy Branch basin in Montgomery 
County, Maryland, The accuracy of the process was verified by 
comparison with aerial photography ground truth. The results 
are summarized in Table 3 , Overall error Was about 6%, ac- 
ceptable for computation of surface cover factors of hydrologic 


models. 


TABLE 3 


Urban 

Forest 

Lakes 

Soil 

Fields 


SURFACE COVER CLASSIFICATION RESULTS, MUDDY BRANCH BASIN 


AREA 

LANDSAT 

% OF AREA A OF 

WATERSHED AERIAL WATERSHED 

INVENTORY 

ERROR 


557 ha 

11 

€49 ha 

12 

-14% 

1.242 ha 

24 

1,292 ha 

23 

- 4% 

30 ha 

1 

29 ha 

1 

+ 3% 

575 ha 

11 

547 ha 

10 

+ 5* 

% it 352 

54 

3,028 ha 

55 

+ 6% 


XV. APPROACH TO THE PRESENT EFFORT 

The current phase of this effort brought the prior work to a 
conclusion through four tasks i 

Task 1 The complete development of the routing model. 

Task 2 The derivation of hydrologic parameters from IhNDSAT 

imagery for an additional watershed, making a total of 
five for which such analysis has been performed. 

Task 3 The documentation of procedures for extracting hydro- 
logic information from IANDSAT imagery. 

Task 4 Demonstration of the procedures and model to a 
selected group of users. 

Figure 1 is a flow chart showing the objectives set out for this 
phase. 

Task 1 constituted the synthesis of a routing module for the 
planning model. This has been completed and programmed in FORTRAN XV 
for operation on small-capacity digital computers . The model features 
user-interactive operation and amenability to present and potential in- 
puts from satellite remote sensors . 

A summary of its derivation is presented in Section VI. (A more 
detailed derivation is included in the Appendix). Task 1 culminated in 
the simulation of a peak runoff event as taken from actual gage records/ 
this is also described in Section VI. 

Task 2 involved the hydrologic land use classification of the 
Bucklodge Branch watershed in Montgomery County, Maryland. The techniques 
developed for Task 3 (described below) were applied to cover 


PLOW CHART OF CURRENT PHASE 








classification of Bucklodge Branch, This task also included the deri- 
vation of model parameters from the remotely sensed information, Re- 
sults are reported in Section VI, 

Task 3 resulted in the development of generalized procedures for 
extraction of quantitative information from LANDS AT imagery through 
simple visual analyses. The findings, given in Section V, showed that 
relatively uncomplicated techniques can yield results adequate for hy- 
drologic modelling. 

The final task included identification of prospective user groups 
and exposure of the model and associated analysis procedures to them, 

A summary of those activities is presented in Section VII. 
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V. PROCEDURES FOR EXTRACTION OF HYDROLOGIC DATA FROM LAND SAT IMAGERY 

The quantities of rainfall which, eventually become stream discharge 
depend in large measure upon the character and condition of the land sur- 
face over which it must run. Rain falling on a watershed is subjected 
to the retardant forces of the surface cover? this retardation, in turn, 
regulates the fractions of the water mass which run overland, infiltrate 
or are evaporated. A rural watershed covered with thick grass 
and trees , for example, will have a much different response to a peak 
rain event than will its urbanized counterpart. 

Figure 2 exemplifies this effect. The curve shows the impact on 
runoff peak of altering the surface complex (as approximated by the S.C.S. 
curve nunfcer) . Points on the curve were generated using a recorded rain 
event input to a 150 square kilometer basin in the region of Tuscany, 
Italy.* The sensitivity is marked. Figure -3, similarly, demonstrates 
the effect on the discharge peak of seasonal change, in this example, 
that of cultivating a portion of the watershed. Again the results are 
pronounced. Two conclusions can be made: 

1) Surface cover is a significant driver of peak discharge 
event, 

2) The provision of periodic land cover information can be 
of value to the hydrologic modeller. 


* This watershed is currently being modeled for the government of 
the Region of Tuscany, Italy. The modeling techniques were derived 
from the present NASA — sponsored effort. 
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Prior work in remote sensing for hydrology has shown that surface 
cover can be derived from presently available LANDSAT imagery. A problem 
with reducing potential to practice, however, has been the absence of 
standard procedures for extraction of required data by uncomplicated , 
inexpensive, yet accurate, means. 

One objective of this effort was the development of procedures for 
utilizing LANDSAT data to determine watershed surface cover. An important 
result is that the procedure eventually developed is applicable not only 
to tiie specific model developed herein, but can be generalized to all 
hydrologic models . 

5.1 Importance of Surface Cover to Hydrologic Models 

For reasons just discussed, all but the simplest hydrologic models 
take surface cover into consideration. Specific implementations vary 
but generally include the following: 

a. Computation of initial detention storage 

b» Calculation of surface resistance to overland flow 

c. Division of the rainfall mass into its overland and intiltrated 
fractions 

d„ Estimation of potential evapotranspiration 

e. Quantification of the effects of impermeable areas 

Figure 4 summarizes the surface cover inputs to several widely-used 
hydrologic models. It is clear that a substantial amount of surface- 
related information is required. 



MODELS 


US DAHL" 70/74 


USGS 

UTAH STATE U. 
STANFORD MODEL 
TEXAS MODEL 
NOAA HYDRO 14 
NOAA HYDRO 17 

API 

SSARA 


C0SSAP.R 

S.C.S, TR-20 
MITCAT 

I 


FIGURE 4 





















For these requirements to be filled by LANDSAT, the frequency and 
resolution of the LANDSAT-derived data must be sufficient to meet model 
demands. With respect to the first point, LANDSAT can provide a new 
surface cover picture every nine to eighteen days . This frequency meets 
and exceeds all the requirements of current models , even when cloud oc- 
cultation is factored in. 

The resolution of land cover detail required by the models varies, 
but is generally consistent with the current capabilities of LANDSAT 
sensors. The exceptions are those models designed to analyze intra-urban 
hydrology? satellite data will be of limited applicability here. For 
those models suited to evaluation of larger watersheds, most require- 
ments can be met through a LEVEL I* land use classification, which can 
be accomplished remotely . The modeller must be able to distinguish 
between: 

a. Urbanized areas 

b. Agricultural (cropped) areas /fallow 

c. Surface Water 

d. Pasture 

e. Forests 


* (As defined in, "A Land-Use Classification System for Use With Remote 
Sensing Data, “ U.S.G.S. Circular 671, Anderson, Hardy, and Roach, 1972) 
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Oils level of detail will permit computation of friction factors 
such as Manning's "n," mensuration of impermeable areas/ calculation of 
surface water area; etc. 

5.2 Alternative Techniques for Analysis of LANDSAT Image ry 

Because LANDS AT data is produced in two forms, computer compatible 
tapes (CCT's) and "photographic" images, two options for its analysis 
exist - computer and visual classification/mensuration. The first 
involves the development of mathematical/statistical algorithms for 
separation of themes based upon the relative output of each of the 
sensors. Visual analysis is the same with the exception that the algorithm 
is mental. The interpreter segregates on the basis of the individual or 
summed reflectance levels of the Multispectral Sensor Bands. 

What are the advantages/disadvantages, if any, of either procedure? 
Our work in hydrology , and other application areas , has made obvious a 
number of factors making visual analysis preferable for "complex" scenes 
and low-cost applications. Principal among these are: 

Training Sample Criteria 

Computer classification and mensuration requires the establishment 
of a number of training samples. These should be "homogeneous" samples, 
containing at least a minimum, number of pixels. With visual analysis, 
it is possible to minimize the training field requirement; it is also 
easier to account for variations within a theme which would result in 
errors with computer classification. 

Maximization of Multi-Temporal Properties of LANDSAT Imagery 

In hydrology, as well as agriculture, the use of multi-date imagery 
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has proved advantageous. Visual analysis facilitates the simultaneous 
evaluation of images from different seasons, years, etc., at low cost. 

Ground Truth 

In most settled areas, the scales of features vary significantly. 

The relatively low, 66 meter average resolution of LANL'SAT causes 
boundary effects which distort the identification and mensuration of 
areas less than several tens of acres, if performed by automatic means. 

Only the visual overlay of aerial imagery or maps against LANDSAT imagery 
corrects this problem. Projection of visual images at scales which 
can be comfortably observed greatly facilitates scale matching. 

Accuracy of Results 

Except for LANDSAT scenes taken at the correct seasons when reflectance 
differences among the various ground covers are most pronounced and in 
which the features to be identified are extensive compared to the IANDSAT 
pixel dimensions, automatic classification exhibits considerable errors. 
Visual analysis employing one or more of the techniques indicated above 
(multi temporal, superposition with aerial photography) , can result in 
significantly improved accuracies . 

Equipment Requirements/Costs 

The initial capital cost of visual analysis equipment (projectors, 
zoom transfer scopes) is considerably lower than that of digital processing 
equipment. For ''complex" scenes, visual analysis has been found in this 
effort to cost on the order of 10% of computer analysis at equal performance. 
It is concluded that at the present state of the art, visual analysis 
techniques are preferable for small users, because they are more likely 
to utilize LANdSAT imagery if the initial costs and skills levels required 


are kept as low as possible. 

Even with simple visual means two problems still constitute impedi- 
ments to the practical use of LANDSAT imagery: 1) the lack of a standard 

procedure; and 2) the unavailability of projection/analysis devices 
amenable to complete analysis of satellite imagery at lew cost. Ecosystems 1 
effort has resulted in techniques for solving both these problems. 

A device for easy analysis was developed and tested* which is suited 
to LANDSAT scales and resolutions. It is described in Appendix B. 

Hie following describes a recommended procedure for deriving 
quantitative hydrologic data from the MSS imagery. 

5.3 Visual Extraction of Hydrologic Land Use from LAI© SAT Imagery 

Hie preceding sections have described the uses of surface cover 
data by hydrologic models, and indicated that many are obtainable from 
LANDSAT imagery. It is significant, in addition, that the extraction 
of these data can be accomplished through simple and inejgoensive visual 
means, procedures similar to those used in analysis of aerial photographs. 
Hie development of a certain amount of expertise on the part of the 
interpreter is necessary, but the level required is not prohibitive. 

Hie most important concern is that a careful and consistent analysis 
procedure be adhered to, such as that described herein. 

5.3.1 Land Use Information Content of the Four LANDSAT Bands 

In order for LANDSAT data to be useful in hydrologic analysis: 

* Hie device was developed on ECOsystems * own funds , but utilizing 
to the fullest the experience garnered in this NASA - sponsored 
effort. 
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resolution of the sensors must be sufficient to provide the required 

* 

level of classification and area measurement} and, the optimm spectral 
bands must be chosen for the surface feature under study* 

With respect to the first criterion, the minimum resolvable area 
of the LANBSAT sensors is of order 0 . 4 hectares (1 acre) * Hie applica- 
bility of LAN DSAT, therefore/ will be optimized where the unit of analysis 
is sufficient, i.e. in medium to large-sized basins predominated by rural 
land use; it is in these basins, though, where gage records are typically 
less detailed and where modelling can have large benefits (snowmelt 
analysis, water works sizing, flood forecasting)* 

With respect to the second criterion, no single LANDSAT band can 

distinguish all land use classes well enough to meet the requirements 

of hydrologic modelling* • However, the output from the four sensors 

taken individually and in sum do provide sufficient information* Each 

sensor sees a distinct spectral range and, therefore, receives different 

radiance levels from a particular object. For the purposes of thematic 

segregation, it is preferable that a sensor see a theme at one extreme 

of reflectance (±»e., very high or very low) while giving the opposite 

value for all other themes. Ibis maximizes contrast and likewise the 

separability of the theme* Again, no single spectral band can accomplish 

this ideal case, but an examination of the imagery shows that hydrologic 

land use class es can be separated. 

■ * 

For example, in LANDSAT Band 7 (0.8 - l.lym) surface water appears 
very dark while its surroundings are generally much lighter whereas 
for Band 5, forested areas appear dark. Figure 5 illustrates this. 
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Table IV summarizes exemplary radiometric data and indicates which 
bands appear preferable for a given theme. Ihe values, derived from 
imagery from the Baltimore- Washington area, show where the land uses fall 
on a fifteen-step gray scale. (A value of 1- 100% reflectance; 15 - 
0% reflectance) . Themes will be most easily separated where the variance 
is small and where the other themes have different values. 

5.3.2 Use of Composite Imagery 

The utility of LANDSAT imagery can be augmented by combining the 
individual bands into composites. The addition of color to each band 
facilitates the separation of surface themes, because it allows the 
analyst to see the information contained in all bands simultaneously. 

A simple, low cost procedure for color encoding individual black 
and white images and compositing them is available in the diazo process. 

Its best us^ge, as developed by ECOsystems, is to produce plastic transparencies 
in shades of a single color directly from LANDSAT 9 n x 9" images. To ac- 
complish this, the black and white image is sandwiched with a sheet 
of diazo film and exposed to ultraviolet light. The film is 


TABLE XV 


Surface Water 

Urban 

Forest 

Fields/Agricultural 


BAWD 4 

OCT. *73 IMAGE 

5-10 
1- 7 
10-11 
4-13 


JUNE ’73 IMAGE 


Surface Water 
Urban. 


Forest 


Fields /Aqricultural 


BAWD 5 

10-14 
1- 4 
. 0-11 
1- 7 


9-10 
1 - 2 
10-11 
1-12 


BAND 6 


Surface Water 

Urban 

Forest 

Fields/Agri cultural 



BAND 7 


Surface Water 

15 

15 

iftfoan i 

10-11 

10-11 

Forest 

7- 8 

2- 5 

Fields/Agriculture 

2-11 

1-11 


indicates best single band for given theme 









then, developed by simple equipment which eposes it to ammonia vapors. 
Figure 6 exemplifies the process product ; a diazo positive of a single- 
band - LANDSAT image processed in cyan. Hie procedure is repeated for 
the other bands producing each band in a different color. Finally, the 
diazo transparencies are overlayed, registered, and glued together to 
form the composite image. Figure 7 is a contact print of the resultant 
composite. 

Seve- 1 advantages to image analysis with diazo transparencies have 
been found: 

1. Hie costs are very low, especially when compared to computer- 
aided classification techniques. Hie exposure and develop- 
ment unit can be purchased for approximately $500 to $1,000. 

Hie materials cost for composite transparencies is less than 
one dollar. 

2. Hie diazo film yields a faithful positive reproduction of the 
original image. Very little detail is lost, and processing 
time is short. A composite transparency can be produced in 
under a hour. 

3. By varying exposure times and film color, images can be 
"tuned" to accent desired features while depressing others. 

Hi is significantly augments the separability of surface 
themes. Several composites can be prepared for a watershed, 
each of which optimizes the visibility of a specific land 
use. 

5.4 Recommended Procedure for the Visual Analysis of LANDSAU Composite 

Imagery 

Hie visual analysis of imagery for hydrologic land use should proceed 



i 



in six steps: 




SECTION OF CONTACT PRINT OF THREE BAND 


DIAZO COMPOSITE, 6 October 1973 


BAND 4 - YELLOW 


BAND 5 - MAGENTA 


BAND 7 - CYAN 
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1) Preparation of "tuned’ 1 Diazo Composites 

2) Image Pro jection/Groun & Truth Overlay 

3) V7atershed Delineation 

4) Analysis of Ground Truth 

5) Segregation of Surface- Coyer Themes 

6} Quantification of Hydrologic Parameters 

5.4.1 Preparation of "Tuned" Diazo Composites 

As noted earlier, no single LANDSAT band contains all the hydrologic 
information; each sensor maximizes the reflectance of the surface in a 
different spectral range. Full advantage of these reflectance differences 
has generally not been exploited in visual analysis; ' rather, past work 
has relied mostly upon LANDS AT 1 s geometric properties - identification 
of shape and form from single color composites. 

The diazo compositing process affords the analyst the opportunity 
to optimize the value of radiometric differences. Consider the fol- 
lowing example: Band 5 (0.6 - 0. 7]iin) receives low reflectance from 

vegetation since these wavelengths are absorbed by chlorophyll. Band 7 
(0.8 - l.iym) , in contrast, receives nearly 50% of the incident radiation 
from mature green vegetation. This phenomenon can be exploited to 
discriminate vegetation. 

Suppose that the watershed under study contains heavily vegetated 
areas. Is it possible to assign a desired color to these areas in 
preparing the composite. 

If it is desired that these vegetated areas appear red, it is only 
necessary to encode the Band 5 diazo positive in red (and/or Band 4, if 
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it is used, since the two are correlated) and to develop Bands 6 and 7 
in contrasting colors. The red vegetated areas of Bands 4 and 5 will show 
through the equivalent sections of Bands 6 and 7 which will he clear. 

“The result is a composite "tuned" to vegetation. 

Thus , by referral to the relative reflectances of land use classi- 
fication, such as shown in Table 4, the analyst can compose his imagery 
to make more obvious the themes he wishes to identify. The technique 
described is valid for themes which provide high contrast (lew or high 
reflectances simultaneously present) . A modification of the process 
allows separation of themes whose reflectances lie closer together (lower 
contrast) . It is a property of diazo imagery that increasing exposure 
time decreases the content of the image; i.e. a long exposure time drops 
out themes wxidi high reflectance (low image density) . A positive image 
can be color-encoded in a given color and exposed to the point where all 
areas with lesser density than those desired are eliminated. This pro- 
cedure is then repeated with the equivalent negative image, but in a 
contrasting color. When the two images are overlayed, the subject then 
will appear in the color produced by combination: surrounding areas 

will be in pure hues. Figure 8 demonstrates the results. 

The two techniques described above separate all but themes with 
very low relative contrast. To further improve the discrimination 
between these, a third option is the use of multi- temporal imagery. 

For example, forested areas and fields can have similar reflectances 
in all the MSS bands. Fields, however, are far less likely to display 
a signature which remains constant over the year. Figure 9 depicts an 
example. Note that in the February image the forest has lost its 





FEBRUARY 1974 

LANDSAT 3 BAND COMPOSITE 


FOREST CONTRAST HIGH 


FIGURE 9 


SEASONAL CHANGE > IN FOREST SIGNATURE 


OCTOBER 1973 

LANDSAT 3 BAND COMPOSITE 


FOREST CONTRAST LOW 






similarity with fields. 

A recommended procedure for integrating these techniques is to first 
develop a standard composite of bands and colors such as Band 4 (yellow) , 
Band 5 (magenta) and Band 7 (cyan) * Hie analyst should then separate 
all identifiable themes; the resultant classification can then be confirmed 
or adjusted by subsequent analysis of composites tuned to a particular 
theme from different seasons. 

5.4.2 Image Projection 

Land use classification to Level 1-2 requires significant magnifica- 
tion of the 1:1,000,000 scale L&NDSAT imagery. This has previously 
been a drawback, of visual analysis hardware - devices which could per- 
mit enlargement to greater than, about 10X have been very expensive. 

She cost can be greatly reduced, however, as exemplified by the imaging 
system developed for such analysis by ECOsystems. Experience with this 
system has shown, that magnifications of about 40 X are convenient for 
land cover analysis, since they permit the direct overlaying of Is 
1,000,000 LANDS AT images upon USGS 1:24,000 topographic maps. To 
insure the color- fastness of the images during analyses, color slides 
should be made from the diazo composite. Hie slides can be made to 
include as large an area as desired, but it is preferable to reproduce 
the diazos at one-to-one scale. The area represented in a 35-millimeter, 
1:1 from a 1:1,000,000 LANDSAT image will be approximately 86,500 hec- 
tares. Larger watersheds can be analyzed in sections. 
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The slides are made with a 35 millimeter single lens reflex camera 
equipped with a bellows. The diazo composite is illuminated from behind 
by a light source matched to the color temperature of the slide film 
used. The photographs are taken with the camera mounted on a copy stand. 
Once exposed, nearly all slide films can be developed "in-house" using 
inexpensive chemistry. The time required from exposure to projector- 
ready slide is less than two hours. 

Finished slides are mounted in the imaging system and projected. 

As shown in Figure 10 , the analysis table is equipped with an overhead 
mirror which reflects the image onto a work surface . The LANDSAT image 
can thereby be superimposed over maps or aerial photographs. A zoom 
control allows the scales of the projected image and ground truth to be 
matched. 

5.4.3 Watershed Delineation /Ground Truth Overlay 

Analysis of hydrologic land cover begins with the demarcation of 
the boundaries of the watershed. Except in regions of very high relief, 
some ancillary information is generally required. This is readily available 
in the form of USGS topographic maps. Stereo aerial photographs, where 
available, can also be used. The watershed ridgs; line is first traced 
on the map. The LA33DSAT image is projected over the map and aligned 
using roads, rivers, lakes, or other prominent features for reference. 

A sheet of tracing paper is placed over both and the ridge line drawn 
on it. The map is then removed. 

5.4.4 Analysis of Ground Truth 

There are a number of widely- available sources of ground truth 
which can be valuable aids to the interpretation of satellite imagery. 



ECOSYSTEMS ULLAGE ANALYSIS DEVICE 


UNIT 


ANALYSIS 
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Among these are: 

1) topographic maps (USGS) 

2) soil surveys (SCS) 

3) aerial photographs 

Topographic maps , though often dated, contain some land cover 
information* Specifically, forested areas, lakes, streams and urban 
areas are called out. Soil Conservation Service maps are most valuable 
to estimation of the subsurface characteristics of a watershed but also 
show streams, roads,, surface water. Many have the soil classes drawn up- 
on aerial photographs which may be used as ground truth if better ones 
are not obtainable. Low altitude aerial photographs are the best source 
of ground truth, but are also the more expensive and least available. 

What are the optimum uses of ground truth , and how would one proceed 
to analyze land cover without it? In the absence of any ancillary data 
whatsoever, it would be necessary to rely totally upon the radiometric 
characteristics of the LANDSAT sensors. Even in this case, it is possible 
for a great deal of land cover information to be culled out. Visible 
themes such as rivers and lakes, some urban areas, and forests can be 
identified by their shape and form. An interpreter with some knowledge 
of the study area and of the spectral properties of land use themes 
could add another level of information. This, in effect, would be the 
same as checking LANDSAT data against a training field. Familiarity 
with the seasonal or phenologic progressions of vegetation is also 
most valuable. 

A competent interpreter can generally adequately classify hydrologic 
land use to Level 1-2 using only LANDSAT imagery. Hie existence of 
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certain types of ground truth facilitates the process and should be 
taken advantage of. Maps and photographs/ even if outdated, can be used 
to verify that an area classified as forest is, indeed, forest. If such 
fact can be confirmed within the watershed, or even in a nearby location, 
the results of these "tests" can be extended over the study area, e.g, 
if a red-blue area in the LANDSAT image is determined to be forest, 
areas of like color can also be classified as forest. The resultant 
classification can be spot-checked where the presence of true forest 
areas is known. 

Prior to evaluation of the LANDSAT images, all available ground 
truth should be collected. The data obtained should be summed into a 
base map for the basin. All surface features discernible from the 
ground truth should be designated on the map. In the ideal case, ground 
truth is sufficiently complete to designate a gross, even though mostly 
outdated, land use map for the entire basin; the function of the LANDSAT 
data, in this case, is to update the older information and to provide 
the important seasonal information which it lacks. 

5.4,5 Classification of LANDSAT Themes 

The next step in surface cover analysis is the actual classification 
of themes from the diazo composite. The image has been projected on the 
work surface and the watershed boundaries established. The analyst must 
then determine how many separable themes, or colors, exist within the scene. 
Figure 11 exemplifies the process. The analyst should make a tracing of the 
themes which are identifiable. Once drawn, it is helpful to assign an ident 
fier to each area delineated and identified* it is suggested that this be 
accomplished through comparison with a set of standard color chips ouch 



FIGURE 11 


SEPERATION OF LAND USE THEMES 


MUDDY BRANCH BASIN 


E2 soil. 
£•3 Urban 



as those used in aerial photographic interpretation- Each separable 


color is thereby given a unique label; this insures repeatability and 
also facilitates validation by other interpreters . 

At this point the analyst has a map divided into color themes 
but has not as yet identified what the colors represent. To do this, he 
refers to the available ground truth; selects a test area from within 
each color class; compares it to the maps or photos , and assigns the 
area to a land use designation as in Figure 12. The resultant land use 
map should be regarded as still tentative. The classifications are 
subject to modification using composites of different composition, i-e- 
different tuning, or from different seasons. Generally, initial classifi 
cation reveals differences between the L&NDSAT imagery and the ground 
truth due to changes more recent than are reflected in the maps or 
photos; these can be confirmed by reference to additional images . 

The result of multiple-image classification will be a hydrologic 
land use map suitable for the derivation of quantitative model parameters 

5*4,6 Quantification of Hydrologic Parameters 

With the land use map in hand, the derivation of model parameters 
can proceed as it would using any source of surface cover data- The 
exact procedure will vary with the model used; the following paragraph 
describes a generalized technique used by several models - the computa- 
tion of an average value of Manning’s "n ir for a subwatershed - 

Planning models often require the division of watersheds into 
hydrologically homogeneous sub-areas . The criteria for defining 


FIGURE 12: CLASSIFICATION OF LAND USE THEMES 
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homogeneity include topography, soil characteristics, and surface cover. 
T/pically, maps of each are prepared and compared in overlay fashion. 

She modeller then divides the watershed to give the optimum uniformity 
among these three parameters. Figure 13 demonstrates the process. 

Once the homogeneous sub-areas are established, an average set of 
parameters is assigned to each. With respect to surface cover, an. 
average surface friction factor, such as Manning’s "n-, " is computed 
for every sub-area. Figure 14 shows how this is accomplished from the 
LANDSAT map as well as from aerial photography. First, the sub-areas 
are delineated on the land cover map. Then the percentage of its area 
covered by each use is computed using a calibrated grid, as shown, or a 
planimeter. Finally, the percentages computed are compared to a table 
of values for Manning’s "n 1 ’ as reproduced in the figure. Their computed 
weighted average becomes the average friction factor for the sub-area. 
Parallel analyses of LAKDSAT and aerial images, performed during this 
effort, yield nearly identical values for surface friction. 

These values become direct inputs to planning models. The following 
section demonstrates the application of the techniques described to an 


actual watershed. 
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VI. DEVELOPMENT OP REMOTE SENSING MODEL AND APPLICATION TO TEST 

WATERSHED 

S.l Model Development 

During the course of this project, several criteria for planning 
models using remote sensing were developed. Ihey include: 

1) The model should consider all hydrologic "drivers, " i.e. , 
tlie surface and sub-surface processes which significantly 
affect the discharge peak rate and duration. 

2) It should take maximum advantage of presently feasible 
remote sensing inputs (surface cover) and be amenable 
to new ones (soils, soil moisture, etc.) as they are 
developed. 

3) A principal function of the study should be to demonstrate . 
to prospective users of LANDSAT data how remote sensing 
techniques can be applied to the general class of plan- 
ning models with actual watershed data. 

A model has been deduced from standard practice and modified to a 
form consis taut with these criteria. It simulates the flow of rain- 
water over the surface of a watershed, subjects it to subsurface abstrac- 
tion. (using the Holtan infiltration formula) ; and routes the excess 
through stream channels to the outlet. Surface cover r land use, soil 
type, soil moisture content, basin, and channel physiography and topo- 
graphy are explicitly considered. 

As presently configured, the model consists of an overland flow 
and a channel flow module. The overland component is based on the 
conservation of rainfall volume: 


“ 4 & 


V =* qt - I 

Where; 3 

V ** overland flow volume r m 

3 . 

q f* rain, rare, m /sec, 

3 

1 * total infiltration volume, m 
t = rain duration, secs, 

The volume of runoff over time is defined by the Manning relation 
for surface flowi 


P * 


h 


tyl 


2/3 


Where » ■ " . ’ 

55 = surface runoff, tn /sec 

9 =» slope, m/m 

n Manning's surface friction coefficient 
y « depth of overland flow, m 


The model operates indiscrete time steps, computing a depth value 
Cyl for the beginning and end of each interval and using their average 
to calculate a discharge (£1. The Q values for all time steps consti- 
tute the overland flow hydrograph, 

The channel flow module takes the overland hydrographs from homogen- 
eous sub-- areas of a Watershed and routes them to the outlet. Manning's 
relationship is used to compute channel outflow. Figure 5 summarizes the 
model's components; a logic flow diagram and a FORTRAN listing are in- 
cluded in the Appendix, Table 5 lists the input data required by the 


model. 

















TABLE 5 

SUMMARY OF ECOSYSTEMS' MODEL INPUTS 


OVERLAND FLOW PARAMETERS 

CHANNEL FLOW PARAMETERS 

Unit strip length, meters 

Channel length, meters 

Unit strip slope, meters /meter 

Channel slope, meters/meters 

Unit strip Manning's "n" 

Channel Manning's "n" 

Initial depth of flew, meters 

Channel Hydraulic radius, meters 

SOIL PARAMETERS 


Soil moisture capacity, meters/meter 
Antecedent moisture content, trie te rs/inete r 
Percolation rate, meters/sec. 

Steady-state infiltration rate, meters/sec 


OVERLAND FLOW INPUTS 

CHANNEL FLOW INPUTS 

Rainfall temporal profile 
Rain rate Cm/sec} 
Time Csecs) 

Unit strip Hydrograghs 

Discharge, or /sec 
Time, secs 


The parameters are input in the conversational mode, i.e,, the pro- 
gram 11 asks" the modeller to enter each .one as required „ the output of 
the overland flow module is printed in tabular form and there is an op- 
tion to graph the data, A. sample run is shown in Figure 16. 


6.2 Application of the Model to Test Watersheds 

Tests were run using LANDS AT data as inputs to verify the operation 
of the -model for typical ungaged watersheds. The tests served two pur- 
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SAMPLE! OUTPUT OF STKCP HTDROGPAPH MODEL 



FIGURE 16 {coat'd} 
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1, To substantiate that the model is capable of producing sat- 
is factory results as compared to similar models of ungaged 
watersheds, The watershed used was gaged and these records 
were used for co;n£>arison, 

2» lb demonstrate for prospective users of LANDS AT data that the 
analysis techniques described in Section IV are applicable 
in ^hands-on" hydrologic modeling and provide sufficient im- 
provement to warrant their use. 

The simulation test encompasses two phases* data collection and 
reduction # and the actual model runs, 

* 

6.2.1 Data Collection and Reduction 

The first step of this phase was the selection of a test watershed. 
After study, the Bucklodge Branch basin, in Montgomery County, Maryland, 
was selected, The advantages of this site were: 

1, Sufficient rainfall and streamflow records were available to 
validate the model's results. 

2, The watershed contains several distinct land uses providing a 
good test of what can be distinguished from LANDS AT, 

3, Aerial photography ground truth existed to validate the re- 
sults of LftNDSAT analysis. 

4, The proximity of Bucklodge Branch permitted any moot points 
to be resolved by field survey, 


54- 


Buck lodge Brandi drains approximately twenty three square kilo- 
meters C9 square miles), of an urban/rural setting in Montgomery County , 
Maryland, The overall stream length is 32 kilometers (20 miles) , The 
predominant land uses are rural and agricultural, including fields, 
forests, bare soil, surface water, and a small urbanized area, These 
land uses are typical of those of ungaqed watershed where satellite 
analysis will have the most value. Buck lodge Branch is monitored by 
the Montgomery County Department of Environmental Protection and is 
equipped with two raingages and a river stage recorder. Figure 17 is a 
map of Bucklodge, showing the locations of the gaging sites. 

The table contained in Section 6.1 listed the data which were re- 
quired before the model could be run. Since the purpose of this pro- 
ject was to rely minimally upon non-satellite data, only data sources 
readily and widely available were called upon, The soils parameters, 
infiltration apd percolation rates, were taken directly from the Soil 
Conservation Service survey of Montgomery County, Figure 18 exempli- 
fies the available soils data, Slopes and lengths were measured from 
a U,S, Geological Survey 1:24,000 scale topographic map. Land cover 
data was derived solely from LANDS AT imagery, The resulting classi- 
fication was verified by comparison with 1:130,000 color infrared 
aerial photography, 

She LANDSAT imagery selected for evaluation was dated 6 October 1973, 
Xt was chosen for image quality and tenporal proximity to the aerial pho- 
tographs CPeceatoer 1973), Multi-spectral sensor bands 4 (0.5-0. 6 jlm) , 

5 (0*6-0, 7 yml and 7 (0, 8-1,1 yml were combined in a three-band diazo 




FIGURE 18 


EXEKPT FROM S.C,S t SOIL SURVEY OF MONTGOMERY CO*, MD* 


Depth to 

seasonally Depth to 
high \vu ter bedrock 
ta 


Brief description of site and soil 


GlenviUe eilt loam, 0 to 3 per- » 1 to 3 
cent slopes* 

Glenvilic silt loam, 3 to S per- * 1 to 3 
cent sloped 

GlenviUe slit loam, 3 to 8 per- 1 1 to 3 
cent slopes, moderately erod- 
ed* 

Gullied litnd, Penn materials {*) 


Huntington »Ut loam, 0 to 3 per- 
cent slopes. 

Huntington silt loam, 3 to 8 per- 
cent slopes, moderately srod- 


Moderately well drained soils in ! 
draws and depressions and up- 
land fiats; developed in residu- 
um weathered from crystalline 
rocks; some colluvial surface 
material. Seepage spots fairly 
common. Fragipan in subsoil. 

Very severely eroded areas that 
originally were Feim silt loam 
or undetermined soils of the 
X 3 etm series; very shallow; in 
many places bedrock is exposed. 

Well-drained, deep soils on flood 
plains; developed in alluvium 
derived from limestone. Infre- 
quently flooded. . 


Depth from 
surface 
(typical 



US IDA 
{textural 
class) 


0 to 5 Silt loam 

5 to 28 Silt Iaain 

28 to 42+ Silt loam. 


0 to 10 ; Chanti cry silt 

loam. 

10 + 


0 to 12 Silt loam 

12 to 00+ Silt loam 


Engineering classification Percentage passing — 


Selected characteristics significant in engineering 


No. No. No* Range in 
Unified AASHO 200 10 4 permeability 

sieve sieve sieve 


Shrink- 

Structure Reaction Dispersion swell 

potential 


ML or CL \ A— 4 or A*-G 


Ifickutt Utfur plf 

0* 33 to 2. 0 Granular to aubangular 4 t 5 to 3* 5 Moderate-^ Low. 

blocky. 

* OG to 2. 0 FJaty and subangular 4, 5 to 5. 0 High Low, 

blacky. 

• 20 to 0. S3 Subangular blacky 4* 5 to 5. 0 High™ Low* 


ML or GMJ A-2 or A-4J 35 50 


GM or GO. 


ML or CL A-4 or A-6 


20 30 


80 100 1G0 

oo lao ioo 


. 20 to 2. 0 Subangular blacky — 4* 5 to 5. 0 High 

* (Fractured shale) 


. G3 to 2, 0 Granular to platy 

1 20 to 0* G3 Blocky and subangular 


0. 1 to G. S 
0. 1 to 7. 4 
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conposite as described in Section 4. The color coding initially se- 
lected was the same as that used by the Earth Resources Observation 
System Data Center, i,e,, Band 4 positive - yellow, Band 5 positive - 
magenta, and Band 7 positive - cyan. Figure 19 is an enlargement of 
the composite, 

The analysis proceeded by transferring the area of the 9" X 9” 
diazo containing Buck lodge Branch onto color slide film. A 35mm camera 
with a slide copying attachment was used; this guarantees the color- 
fastness of the image to be analyzed. The slide was projected using 
the Mark I Ecosystems Universal Image Analyzer at a 1; 24,000 scale, 
(.41,7 X magnification!. 

The specific steps required to translate the colors to specific 
land use and coyer classes were described in Section V and are sum- 
marized in the following! 

Step 1; Location of Watershed ; The ridge line of the Bucklodge 
Branch watershed was first demarked on the U.S.G.S. topographic 
map. This outline defined the area to be classified. The Eco- 
systems analyzer permits magnification of a projected image to 
predetermined scale. By adjusting the projection distance, the 
LANDS AT composite was enlarged to match the topo-map scale. The 
image is reflected off an overhead mirror and superimposed on 
the base map, An exact superposition is attained by adjusting 
the zoom lens until alignment is achieved between lakes, roads 
and/or other prominent unambiguous features in the scene. The 
watershed is thus bounded by the ridge line drawn on the map. 


FIGURE 19: LANDSAT DIAZO COMPOSITE, BUCKLODGE BRANCH 
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A sheet of plain paper was placed over the topo map and the basin 
boundary traced on it, The map was no longer needed and was re- 
moved* 

Step 2 : Preliminary Land Use | Classification t An initial assess- 

ment of how many separable themes exist within the basin was made 
by visual interpretation, ’’Themes" in this sense actually refers 
to a color in the LANDSAT composite. The interpreter determined 
how many separable colors could be identified by reference to a 
set of standard color chips as indicated in Section V* 

Step 3 1 analysis of Ground Truth » Aerial photographic ground 
truth was used to test classification accuracy of the LANDSAT 
imagery. It was prepared in color slide format and projected, in 
fashion similar to the savellite composite, to-a scale of 1i24,00Q. 
The watershed ridge line was marked off. In this case, however, 
the interpreter identified individual land use classes rather than 
color. The resolution of the aerial photography was sufficient to 
permit separation to Levels 1-2 which is adequate detail for most 
modeling applications i Five classes were separated - forests, 
fields ,1 agricultural Cf allow) f surface water, and urban areas. 

Step 4 1 Correlation of Composite Color and Land Use Class ; Exam- 
ination of selected test areas in the LANDSAT image with corres- 
ponding points on the ground truth photograph yielded a correla- 
tion matrix of color and land usei 


COLOR 


LAND USE 


Black 

Blue /White 
Light Red 
Dark Red 
Red/Blue 


Surface water 

Agricultural (Fallow) 

Fields 

Forests 

Urban 


Step 5 1 Analysis of Classification Results : Hie area of each 

class was measured with a calibrated grid overlay. The results 
are summarized in Table 6. 


TABLE 6 

SUMMARY OF CLASSIFICATION RESULTS 


CLASS 


Fields 

Forests 

Urban 

Surface Water 
Agricultural (Fallow) 
TOTAL 


AREA, LANDS AT AREA, AERIAL 

IMAGE (ha) PHOTO (ha) 


RELATIVE 

ERROR 


79 

1,510 

-2 

61 

738 

+ 4 

10 

9 

+n 

7 

8 

-17 

87 

85 

+2. 

50 

2,350 

2,7 


The classification results indicate a high overall accuracy. It 
was observed that large (greater than 10%) relative errors resulted 
only from classes with very small absolute areas (less than 1% of total 
watershed area), In the case of urban areas, for example, the 9 hec- 
tare total area was distributed over six separate sites, none larger 
than about three hectares. The total surface water area of 3 hec- 
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•tares was divided among nine individual lakes. Figure 20 indicates that 
the relative error is a monotonically decreasing function of area of 
the class, as supported by theory. All points fall close to the theo- 
retical line with the exception of agricultural (fallow) areas, which 
displayed an error less than that expected, This further indicates 
that, hecause of the radiometric characteristics of agriculture lands, 
October is a good time to "look" for these, i,e, , their contrast with 
surrounding areas is particularly high, This demonstrates the value 
of multi-temporal imagery indicated in Section V . 

6.2,2 Peak Event Modeling of Bucklodge Branch 

A simulation for all of Bucklodge Branch was not attempted since 
the stream gage is not located at the outlet, Bucklodge Branch was 
modeled down to the gage point in order that the results could be veri- 
fied, Figure 21 is a reproduction of the U.S.G.S. 1:24,000 scale topo- 
graphic map of this section. The heavy lines delineate inport ant: sub- 
basins, The data collected earlier indicated that there were no dra- 
matic variations in soils, topography, flow length, or surface cover 
within any subbasin. 

The model used approximates the natural basin subareas with rec- 
tangular strips possessing averaged surface and subsurface character- 
istics, For Bucklodge Branch, each subbasin was replaced by two strips 
representing the areas on. either side of the main channel as shown in 
Figure 22. The area of each strip was set equrl to that of the subarea 
it represents. Its length is assumed equal to the channel length. 

This fixes the width, which is indicative of the overland flow length. 



ABSOLUTE ERROR (%) 


FIGURE 20 


LANDSAT vs AERIAL NVENTORY ERROR - BUCKLODGE BRANCH 










Slopes and soil parameters are averaged over the area of the strip. 
Table 7 summarizes the parameter assignment for each strip. 

Figure 23 presents the hydrologic schematic for the surrogate 
watershed, She modelling process now requires that rainfall from a 
test event be input to each strip and the resultant runoffs routed 
down the channels. Characteristic rainfall test event were selected 
from among records obtained from the two rain gages in the watershed, 
The rain and resultant runoff stage readings are reproduced in Fi- 
gure 24. .Rainfall inputs applied to each strip were taken as a 
weighted average of the two raingages. 

Model runs for the particular rain event selected showed that no 
overland flow was produced by thirteen of the twenty strips. Over- 
land flow from the remaining seven were input to the routing module 
and the outflow hydrograph was derived* Preliminary runs indicate 
that overall basin response is heavily dependent upon the initial 
soil moisture condition assumed; a reasonable level of initial soil 
moisture was set, therefore, based on the analysis of rainfall for 
14 days prior to the event* It should be noted that sensitivity to 
antecedent moisture is an important common characteristic of planning 
models r it may be a source of additional benefits from remote sen- 
sing of soil moisture in the future* 

Results of the overland and routing runs are shown in Figure 25. 
She model exhibited a good approximation of actual discharge as corn-' 
puted from the gage records and rating curve* The results to date 
are encouraging* Even more important from the standpoint of modeling 

















TABLE 7 


BUCKEODGE - H7DR0L0GIC DATA OP COMPONENT STRIPS 


No. 

STRIP 

LENGTH 

(M) 

WIDTH 

(M) 

AREA 

(M 2 ) 

SLOPE 

MANNING’S 

N 

MOISTURE 

CAPACITY 

(H) 

PEHC. 

(M/SEC) 

FINAL INFILTRATION 
(M/SEC) 

I 

AIL 

402.34 

1626 

6504.00 

.040 

.076 

.038 

. 00000134 

.000004445 

2 

AIR 

271.27 

1192 

323353.84 

.060 

.076 

.036 

.00000119 

.000004445 

3 

AR 

463.29 

1051 

486917.79 

.060 

.044 

- .037 

. 00000105 

.000004445 

4 

AZL 

509.02 

847 

431292.40 

.060 

.031 

.039 

.000001129 

.000004445 

5 

BL 

505,97 

852 

431086.4 

.070 

.045 

.063 

.000001917 

. 000004445 

6 

BR 

298.70 

1207 

360530.90 

.080 

.045 

.036 

.00000105 

.000004445 

7 

CL 

384.05 

1045 

401332.25 

.080 

.045 

.035 

.00000184 

.000004445 

8 

CR 

408.43 

1237 

505227. 90 

.080 

.057 

.038 

. 000001129 

. 000004445 

9 

DL 

542.54 

719 

390086.26 

.050 

.056 

.039 

. 00000127 

.000004445 

10 

DR 

304.80 

793 

241706.4 

.060 

.037 

.060 

.000001199 

.000004445 

11 

EL 

353.57 

841 

297661.70 

.060 

.055 

.056 

. 000001129 

,000004445 

12 

ER 

323.09 

518 

167360.62 

.060 

.062 

.066 

.000001341 

. 000004445 

13 

FL 

548.64 

955 

523951,2 

.060 

.044 

.052 

.000001199 

. 000004445 

14 

FR 

201.70 

443 

89353.10 

.090 

.045 

.030 

.000000847 

.000004445 

15 

GL 

262.13 

425 

111405.25 

.130 

.084 

| .070 

.0000014110 

. 000004445 

16 

GR 

457,20 

528 

241401,60 

.120 

.081 

.065 

. 00000134 

. 000004445 

17 

HL 

286.51 

986 

282498,86 

.110 

.084 

.066 

.00000134 

.000004445 

18 

HR 

310.90 

1100 

341990.00 

: .070 

.086 

.048 

. 00000119 

. 000004445 

19 

IL 

329.18 

1230 : 

404891,4 

.110 

.075 

.066 

.00000134 

. 000004445 

20 

IR 

280.42 

967 j 

j 271166,14 

.110 

.097 

.047 

.00000134 

.000004445 



FIGURE 24 i BUCKLODGE BRANCH WATERSHED PRECIPITATION AND DISCHARGE 
GAGE RECORDS EVENT OF 15-17 JUNE 1976 
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FIGURE 25 


KUKF ALL/RUNOFF EVENT OF 12 - 17 JUNE 1974 
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however, is that the analysis techniques developed have been success- 
fully related to actual hydrologic situations, This applicability 
extends beyond the model developed in this effort , and should prove 
valuable to a large class of modelers and users within the hydrologic 
community, 

6.3 Implications of the Introduction of Remote Sensing Data into 

Hydrologic Planning Model 

1, Satellite remote sensing can yield assessments of hydrologic 
land use sufficient for modeling purposes. The radiometric 
and multi temporal properties of the imagery can facilitate 
classification and provide new and additional information 
with respect to seasonal and developmental changes, 

2, This information will be of value to a large section of hy- 
drologic applications, 

3, The assessment of soils and soil moisture conditions is pre- 
sently difficult for modellers due to its temporal variabil- 
ity and the lack of records. Future inputs from remote sen- 
sing in this area would be most valuable, 

4, The techniques derived and the demonstration made herein 
should function as a guide for the application of LMDSAT 
data by hydrologic modellers to their models on basins of 
concern to them- 
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VII. TECHNOLOGY TRANSFER TO THE USER COMMUNITY 

The preceding analysis has indicated the applicability of remote 
sensing from LANDSAT to hydrology by confirming its adequacy as a 
source of surface cover information t The remaining problem was how 
best to introduce these results to users, Though hydrologists are 
generally familiar with conventional (aerial) remote sensing, the 
satellite data stream has not as yet seen wide use. This section will 
describe the problems involved in the technological transfer and the 
results achieved to date, and will present recommendations for future 
efforts. 

7,1 Description of the User Community 

Results from a previous NASA-GSFC/ECOsys terns effort, synthesized 
following, characterized the hydrologic model users: 

1, In the United States , agencies and organisations concerned 
with water resources number in excess of 6,000. 

2, The principal contributor to research, development and imple- 
mentation of water-resource models is the federal Government, 
with over 80% of the budget, Most of the federal effort is 
carried on by 11 federal agencies and their contractors, 

3, In descending order of activity are state agencies (several 
hundreds) f state water resources institutes (50) , universi- 
ties (.70 principal) f local governments (in excess of 3,000) , and 
private contractors {.approximately 3,000), 
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4, State water resource institutes and universities are pri- 
marily funded by the federal government. Mast private 
contractors primarily support local Governments. 

5, The eleven federal agencies which primarily "drive" hydrologic 
efforts are | 

Department of Commerce - National Oceanographic a 

Atmospheric Administration 

Department of Agriculture 
• a. Agricultural Research Service 

b. Soil Conservation Service 

c, forest Service 

Department of the Interior 

a. Geological Survey 

b. Bureau of Reclamation 

c. Pish and Wildlife Service 

d. Bonneville Power Administration 

Environmental Protection Agency 

Department of Defense ~ Army Corps of Engineers 

Tennessee Valley Authority 

6, The hydrologic community is diverse in its membership and in 
the models they use. There exists, therefore, no single focus 
at which satellite data analysis can be integrated into stan- 
dard practice. 
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7, The diversity of users has recently been expanded with the 
emergence of several federal water resources planning pro- 
grams which tend to move some of the activity away from 
federal groups and toward local and regional agencies, 

This further accentuates the need for using single analysis 
procedures and hardware that can be employed within the con- 
straints of local government budgets, 

7,2 Impediments to Expanded Use of Remote Sensing in Hydrology and 

Suggestions for Overcoming Them 

In addition to the diffuse make-up of the user community, two other 
problems are experienced in amplifying the use of LANDS AT data. There 
exists in any research activity a natural applications delay between in- 
vention and acceptance. New procedures must compete with everyday prac- 
tice and must overcome an understandable inertia. In the case of 
LANDS AT data there exists a certain skepticism not dissimilar to the one 
encountered by any innovation, it is increased by the limitation of re- 
sources available to any one user, which makes him naturally hesitant to 
attempt new techniques without a certainty of their results, Such skep- 
ticism is healthy, however, and forces that the claims made for L&NDSAT 
be thoroughly tested and made relevant to the user, 

A recommended strategy for exposure of the new procedure, therefore, 
is not to represent LANDSAT as a panacea for all problems? rather, speci- 
fic applications of interest to specific users should be "benchmarked, " 
For this reason, this effort did engage in a specific benchmark test, 
as reported in Sections V and VI of this report, The test and the pro- 


cedures developed and described herein can be repeated by any small or 
medium user on his own watershed, and the accuracy and cost of results 
thereby certified for his specific application! 

7.3 Summary of User Contacts Made to Date 

Efforts made during this study to foster the expanded use of re- 
mote sensing for hydrologic modeling encountered the impediments des- 
cribed above. Nevertheless, significant successes were also encoun- 
tered with overseas users. Principal among these werei 

1, A hydrologic modeling effort, based essentially upon the de- 
velopments accomplished under NASA sponsorship, was contracted 
by the Region of Tuscany, Italy* 

She model embraces a watershed of approximately 90,000 Hec- 
tares, with mixed urban, agricultural, forest, mountainous 
contents, The model must be able tos 

al Predict runoff events in real tine for matching supply 
against demand. 

bl provide flood warning 

c) provide accurate runoff predictions for unusual events 
(planning model) for the sizing of waterworks 

dl Simulate the effect of changes in land use — reforesta- 
tions, industrial parks, urban expansion— upon the 
watershed's hydrologic regime. 

She effort additionally includes complete transfer of 
modeling and remote sensing technology to the customer. 
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This effort proved that the technique of planning models 
developed herein could be extended to general hydrologic 
modeling. The Tuscany model includes quarterly updates of 
the ground cover of the entire watershed from LANDSAT-de- 
rived information, 

(■> 

2, The specification of a facility for the analysis of LANDSAT 
data was contracted by the same customer. This includes vi“ 
sual and computerised data processing, employing their on- 
premise 370/135 computer, and using the LANDS AT data derived 
from the Telespazio receiving station located near Rome, 
Italy, The LANDSAT data are to be used for general land-use 
information as well as for hydrologic modeling. 

3, An image- analysis device, of the type described previously, 

/ 

was acquired by the University of Naples, Italy, for the 
analysis of land use, water pollution and the inventory of 
selected crops. 

4, A land-use map of high precision, including hydrologic pro- 
perties, was contracted by the Region of Basilicata, Italy. 
The techniques employed will be those developed under this 
effort and reported herein. 

The high interest of the foreign user community is attributed to 
the relatively small number of agencies working in hydrologic modeling 
and the concentration of these activities within regional governmental 


users , 
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The domestic situation differs in that the level of the techno- 
logy transfer effort required is proportional to the number of users 
involved. On the domestic scene, initial contacts have been made 
to date with U.S ,G«S v , Water Resources Division; US DA, Agricultural 
Research Service, and the Maryland Department of State Planning. As 
previously discussed, cultivation of interest with these and other 
organizations will he greatly aided by concise benchmark tests. 
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TABLE AX 


FORTRAN LISTING OF STRIP HFDEOGRSPB MODEL 


LENGTH, 132 ... . . . 

DIMENSION 9t53riiRl5g«h3!j«).F«50g)rXLHU3«hLKZ5hV(E5).Krz5)»CrRU5} 
ra 10 FORMAT CENTER LENGTH (M). SLOPE* MANNINGS H') 

133 n FORMAT (F7.2.2F5.4) 

SO 33 FORMAT CENTER TIME STEP INCREMENT (SECS! AND NO. OF TINE STEPS'! 

1-0 3* FORMAT CENTER TIME SCALE STEP (SECS)’) 

.-.3 W FORMAT (F4.8iI3) 

.78 i5 FORMAT ('ENTER TOTAL *I«E lIMIT (SECS}') 

46 FORMAT (16! 

■31 94 FORMAT UilthZh' TIHE(SEGS)M3X. * S ISCHARCE 4H3/SEC* I S -4)', 

:3Z ^ pqrjsAT (/*7X*l0tF5.2(5X)iF5iZ) 

IS? Z05 FORMAT (5X.F?.Z,5X t Fi5,13*5X,Fi5 l ig*5XiFL5.t0*5X ! Fl5.l0i5X.F15.lO) 

193 53 FORMAT CENTER NUMBER OF RAIN INCREMENTS' ? 

230 60 FORMAT (13! 

213 7 S FORMAT CENTER RAIN RATE (M/SEC! AND DURATION (SECS) *1 

128 m FORMAT C TOTAL RUN TIME MUST BE ACCOUNTED FOR') 

ITS 73 FORMAT (FB.7tF5.91 

1! FORMAT (2> c 6.4*2Flt.9) 

353 72 FORMAT C5X* • * 3X t » I * 1 180A11 

2'5 99 FORMAT (///.38XC+ INDICATES INFILTRAT1O?)',/*30X.'» INDICATES DISCHARGE'./ *33 XC# INDICATES RAINFALL'! 

Itn FORMAT CENTER SOIL MOISTURE CAPACintHli ANTECEDENT MQIS-JREiH) * PERC(«/SEC!i FINAL U1 FIl.(H/SEC!»’. 

122 14 1 FORHAT (3F9.7) 

230 15i FORMAT t /7 / // »6X» »TIHE4SECS) * »7X i ' RfllH(H3/SEC} 1 t7Xt r INFIL. CH3/SEC) T 1 US. 'HtTil * t i3X t *FL0UiH3/S£C) * * 1SX t 'GAP- tHl * ) 

290 HSITE (108*13) 

300 READ (135*20) XL.SL.EN 

US 0RITE (108*30) 

;?8 READ (105*40) DliU 

..J WRITE (108,45) 

?4.< READ (105*46) N! 

,50 NRITE (108*50) 

READ (105,60) 01 
•:§ WRITE (138*181) 

'M READ (105*91) CAP.APtfPERC.F3 


i;i 

M(i)=0 


430 

WRITE 1108.70) 
WRITE (180.80) 


4?(J 

DO 102 I-1.J1 


* * * 

’«♦ READ (185.93) V(I)iJU) 


a.} 

J2=.il+1 


450 

00 til 1=2*J2 


*60 

ill «(I)*IDINT(UC-t)/Dl)+rt(Mi 


478 

BO 121 L“2+ J2 


463 

LI : 5(1-1) + 1 


490 

a * mu 


180 

DO 121. K1=L1*l2 


510 

12! R(K1) = VU-il 


523 

Ml - HU2) 


5^ 

N S IDIN 7 (N1/Dl) 


- , 1 

F5=(SL«.3/EN) 


950 

WRITE (108,151) 


560 

0=3. 


570 

SDELT=0, 

ORIGINAL 
OF POOS- 

588 

590 

K=0 

id) * CAP-API 

595 

F3-F3H41732. 

683 

DO 321 1*1 rN 



TABIiE AI (cont T d) 


K-K+l 
• 8 Di< I 

y (H ,LE, 0.1 GO TO 251 
(R(I>.LE,0.) GOTO 231 
231 IF (S(l) .Li. 0.) GO TO 271 
Sdl - $d)*39.37 
F < 11 = t(StI)**l.*+F3)l 
Fill 8 F?!)/ 141732, 

SU) - S (I*/ 39.37 

IF (H .LE. 0.) GO TO 281 
241 8(1) - F5*.187*H**1,67 
GO TO 291 

2:1 IF (3(11 .LE, 0.1 GO TO 255 
IF (3(11 .LE. 0.1 GO TO 257 
$(Z) 8 Slll*39.37 
?<I) = ((S'I)**l. 4*F3)1 
StU 8 3(11/39.37 
fill 8 FtP/ 141732. 
if (F(D .GT. 8(11) Ftli 8 Rtl) 

253 IF (K .LE. 0.1 Gp TO 281 
QU) = F5*.307*H«U7 
GO TO 29 i 
255 F(n =g. 

P 8 0. 

Sd) »f. 

GO TO 291 
257 3(11 8 0. 

Fill 8 P ERC 

IF (Fill .GT. Rim FCIi 8 Rtl) 

H 8 t.o*Bl*(R(l)-F(I)> 

GO TO ?53 
27 l 6(1) s 6. 

Fdl 8 PERC 

ri 8 H+(1.6*0t+(R(I)-(Q(I-l)/XU“F(I): 1 
GO TO 253 
281 H 8 0, 

Qd) s 0. 

291 SDELT 8 (F(l)-PERC)*B1 
IF (30E:T .LT. (StIKAPJl SDELT 8 Sdl-OftP 
Sn+t)-S(Il -SDELT 
•F (Stmi.LT.g.l 3(1+11=0. 

IF (SU+il.GT.CAP) 3(Id)=CftP 

IF (K .LT. J) GO TO 321 

WRITE ( 188 i 205) TfRtll .F(II iHattCIl iS(Z) 

K=0 

321 CONTINUE 

DATA ST AR . BLANK r DOT . S,_SHt PLUS . PND/ ' i 1 V.’.dV+V? 1 / 

37 FORMAT (2X.F5.0.3I.A1 . 190A11 
WRITE (108.341 

34 FORMAT ./////.’ESTER MAXIMUM DISCHARGE - IX.XKM 
READ (185.38! D 

38 FORMA" (FU.101 
WRITE (198.39) 

READ (105.60) 12 
II 8 I2/D1 
WRITE (108.99) 

WRITE (108.941 
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% 

1 

H 


| 

j 

1 

i 

1 

l 

t 


j 

J 

f 

i 

t 

i 

I 


83- 


TfiBLB AX (eont'd) 


1H# Hi=0jrt!Jg. 

1150 Ot=8 

1163 CTR(l) -- 0, 

1200 TIKE =3, 

1210 D'2=H1*18. 

1220 DO 301 Mi 10 

1236 301 CTRII+il = 1*02 

1232 WRITE (103) 95) tCTR(Ilrt=!>ll) 

1235 00 308 J=lii00 

1237 308 XLN(J) = DOT 

1250 00 309 1=10.100)10 

1260 339 TIN (I) = 3L3H 

1270 WRITE (108i92) XLN 

1236 DO 501 l=lXtN,U 

1290 fiKE = TIKE+I2 

1295 Fill = F(U*XL 

1298 R(I) = RUML 

1299 S=$LSH 

1300 .jAdOIKK (F(l) *101801 /HI-*. 5) 

1301 iF U6 .EQ» 0) G=PLUS , 

1304 05 = IDINT{(R(il*lS030,)/Hl * .5! 

1305 IF U5 .ES. 3) G=PRD 

1308 J3 = IDINTUG(l)+18000.)/tU +.5) 

1309 IF (58 .EG. 0) G'STftR 

•310 LN=lH+1 

1320 IF (LK .EG. 51 GO TO 315 

1330 BO 103 JK= 1 i 188 

1340 103 XUi(JK) = BUNK 

1350 104 DO 105 JK=10. I80il0 

1368 105 XINUK1 = SLSH 

1363 IF (J6 .EG. S) GO TO 504 

1365 XUU-J6) = PLUS 

1380 504 IF (55 .EG. 0) GO TO 508 

1385 XLIUJ5) = PND 

1386 503 IF (J8 .EO. 01 GO TO 510 

1387 XUUJ8) = STAR 

1398 ‘ 510 WRITE (103,37) TIKEiGiXLN 

1400 GO TO 501 

1410 315 DO 189 JK=lit93 

1428 109 XLNUK1 = DOT 

1430 LH = 0 

1440 GO TO 184 

1458 501 CONTINUE 

1466 STOP 

1470 END 
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MULTI SPECTRAL IMAGE ANALYZER 


SYSTEM FEATURES 

The system is designed to perform scientific analysis of imagery j 
in particular from Landsat and Aircraft. It embodies in a single im- 
plementation the features found in existing image-analysis equipment, 
plus additional capabilities, The system is composed of an Image Pro- 
jection and Control Unit which includes t 

1) The Image Former 

2) The Image Work Table unit, 

The system operates in two basic modes i a) tire Image Combin- 
ation Mode f bl the Image Compositing Mode. 

In the Image Combination Mode, up to four black and white trans- 
parencies of the same scene are inserted into the Image Former, The 
four images are combined optically and projected simultaneously -onto 
the Work Table as a single registered image. Various operations can 
be performed independently on each of the Cup to four) images: sel- 

ective filtering tcolorl , flickering, fade-in- fade- out, light inten- 
sity dimming. The scale of the combined image can be magnified up to 
100 times, allowing for example the superinposition of satellite ima- 
gery onto large-scale aerial photographs, topo maps and similar. 

In the Image Compositing Mode , the black and white image repre- 
senting each band is first developed onto a monochromatic image, which 
reproduces the original to a high fidelity of geometric registration. 
This can be done in any one of 9 color-hues. It is then possible to 
operate in two modes j 


a) Inserting each monocolor transparency (up to four) , repre- 
senting a different band, into the Image Former, superim- 
posing and registering them, and projecting the combina- 
tion onto the IVorA: Table, 

b) Superimposing by simple means (external to the system) two 
or more single color images, representing different bands 
of the same scene, and constructing therefrom a single 
multicolored image. This image can then be inserted in any 
one of the four Image Former projectors and projected onto 
the Ft fork Table for analysis. 

Method b) allows the user to construct his own color composites, 
up to 220mm format, with the advantage that he can vary the density 
and color of each image to bring out the most desirable features. The 
economic advantage is the very low cost of compositing, which allows 
ample leeway for experimentation and optimization. The quality is 
fully comparable with EROS-type products, 

A special case of working with composite images is the use of al- 
ready-prepared color transparencies such as those purchasable from the 
EROS Data Center, 

The basic advantages of the system are the followings 

It can precision-expand the images up to 50:1, and up to 100:1 
at lesser precision, thereby allowing superimposition of a wide 
range of imagery scales of the same area. 

It provides the analyst with a clear and large image on a con- 
venient, solid work table, without forcing him to look through 
lenses or optical devices, thereby providing the best match 
between man and machine. 

It provides rapid, smooth, remote and flexible capability for 
registering images, changing colors, and applying various pho- 
tointerpretation devices, allowing the analyst to use simul- 
taneously all the tools available for maximal extraction of in- 
formation from imagery, 


It combines into a single device the features of the two basic 
systems currently available on the market, i.e, » 1) the trans- 

fer type systems, which project an imaqe onto a flat surface, 
allowing superirrroosition ; and 2) the color-additive systems, 
which combine various images of the same scene onto a single 
image. Existing transfer systems within the same price range 
do project single images , but do not combine multiple images t 
existing color- additive systems combine images, but have li- 
mited projection capability. 

The system provides the added capability of reproducing high 
quality, density-controlled, precision-registered color im- 
ages from the original black and white, at very low cost (ap- 
proximately $0,25 per copy), This capability represents an 
added degree of freedom in the radiometric domain. 

The salient items of comparison between the systems are shown in 

Table Bl. 

hn important consideration is the size of the image acceptable by 
the system, Landsat and most aerial images fall prevalently into the 
sizes of 70mm C54mm information size) and 9” X 9" or 225mmm (185mm in- 
formation size). In the case of Landsat, both sizes contain the same 
information t except that the Landsat pixel size is approximately 24 
microns in the 70mm, 70 micron in the 220mm format. Since the eye's 
resolution is of order 70 microns , the 70mm imagery must be expanded 
5-10X and the 220mm imagery 2-3X to reduce the contribution of the eye 
to within acceptable limits. 

The expansion can be accomplished either by microscopy or by pro- 
jection, Experience has shown that projection matches far better the 
human analyzer's ohysiopsychological preference. This is the first re- 
quirement for magnification, The second requirement stems from the 


need to match the Lands at scale with commonly available ground truth 

scales. This is where the high magnification is most important. 

Experience has shown that very seldom does one need to work with 

2 

the entire Lands at frame £35,000 Km )« The overwhelming majority of 

applications deal with detailed studies of areas of hundreds, at most 

2 

a few thousands of Km . Ihe use of the entire Lands at frame is mostly 
limited to sele cting the smaller area to be studied. Hence, there 
exists very little practical need, and significant added cost, to pro- 
vide the capability to magnify 50 or 10 OX the large 220mm image. In- 
stead, the Multispectral Analyzer System allows an equivalent capa- 
bility to the user*, namely to copy the 220mm image in essentially 

perfect rendition, and to cut therefrom the selected 70mm area (re- 

2 

presenting 3,000 km 1 for projection* 
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o Handles up to four image transparencies of format 70mm., si- 
multaneously < 

o Each, of the four images can be spectrally filtered indepen- 
dently, via insertion of any combination of three quality 
color filters (Wratten #25 Red, #47 Blue, #8 Yellow, or 
equivalent! per image, Open gate position also provided. 
Filter insertion independently and remotely controlled on 
all four images (maximum of thirty two different color com- 
binations! , 

o Registration of all four images, or subsets thereof, by se- 
parate, remotely controlled three-degree of freedom adjust- 
ments; rotation, up-down (Elevation), right-left (Azimuth), 

Caution; the capability to register the four images onto 
each other depends upon the precision of manufacture of the 
original imagery. For example, Landsat four-band imagery 
of the same scene is generally registered adequately, Landsat 
images of the same scene taken at different times (multitem- 
poral) are in some cases not exactly registered one with the 
other, In this latter case, the equipment permits selective 
registration of any portion of the image under study, 

o Brightness of each of the four images vary independently 
by means of remotely controlled rheostat type dimmers. 

The brightness- range is continuous from zero to maximum il- 
lumination, This permits the maximum separation of color un- 
hampered by contrast/brightness effects, 

o Variable flicker capability remotely controlled with full in- 
terruption of light on all four images. Two flicker modes are 
provided; manual remote and automatic remote, In the auto- 
matic mode, the flicker continues automatically until turned 
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off by the operator: the flicker frequency is variable 

remotely by the operator. In the manual mode, the opera- 
tor selects one or more of the images: flicker is pro- 

vided by depressing a foot pedal. Typical useful flicker 
frequencies empirically determined are of order 2 cps . 

o Independent and remote focus on all four images. 

o Overall scale expansion on all four images, ratio up to 

1:50 {1:1,000,000 scale imagery can be expanded to 1:20,000) 
in precision mode, up to 100:1 maximum. The scale range is 
set manually by translation of the image-forming equipment. 
The fine adjustment is provided by a remote-controlled, ap- 
proximately 1:1.5, continuous scale expansion on all four 
images . 

o Manual fade-in fade-out control of the image by means of 
separate control which varies the illumination ratio be- 
tween the imagery under study and the superimposed aerial 
photography or map. 

o Capability to photograph the image projected onto the work 
table. Tliis is provided by means of a sturdy camera mount* 
easily attachable and detachable. 

o A vacuum system under on-off control by the operator can 
assist in holding the tracing paper or maps firmly onto 
the work table. 

The control panel is of the "remove-and-relocate" type. 

It is movable to several positions (e.g., on either side, 
or one left one right of the operator) , to best comply 
with, the operator^ preferences. Mechanical relocation 
of the control panel is accomplished by means of pre- 
arranged brackets and corresponding fasteners*, electri- 
cal connections are relocated by means of multiple plugs. 

o Maximum dimension of image projected on work surface: ap- 

proximately 60 x 50 centimeters. 



o Precision quality color-corrected lens train (Buhl variable 
focal length or equivalent) on each projecting system- 
Registration accuracy better than 20 microns at image center, 
no worse than 40 micron at edges (Referenced to original 
image) * 

o The capability to evaluate image errors and distortions is 
provided by means of separately inser table optical reticle 
bearing calibration patterns* 

o Ruggedized metal construction, Operating ambient temperature: 

14°C to 27°C. Storage temperature: -40°C to +50°C. Shock 

and Vibration: Laboratory environment- Aero or auto transportable 

when packed in box supplied with equipment- 

i 

o Pcwer input: 110 volt A.C., 50-60Hz, angle phase approximately 

2,000 watts. 

o Fuzing: Reseattable power breakers. 

o The equipment is designed for easy access for maintenance 
and repair without the need of specialized tools. 
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